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SUMMARY 


A set of nonlinear equations of equilibrium for an elastic wind 
turbine or helicopter blade are presented. These equations are derived 
for the case of small strains and moderate rotations (slopes). The deriva- 
tion includes several assumptions which are carefully stated. For the 
convenience of potential users the equations are developed with respect to 
two different systems of coordinates , the undeformed and the deformed co- 
ordinates of the blade. Furthermore, the loads acting on the blade are 
given in a general form so as to make them suitable for a variety of 
applications. The equations obtained in the present study are compared 
with those obtained in previous studies. Finally, it should be noted that 
this report represents the first in a series of three reports documenting 
the research performed under the grant. The second report (UCLA -ENG-7880) 
deals with the aeroelastic stability and response of an isolated horizontal 
axis wind turbine blade. The third report (UCLA -ENG-7881) deals with the 
aeroelastic stability and response of the complete coupled rotor/tower 
system simulating essentially the dynamics of the NASaA)OE Mod-0 config- 
uration. 
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1. INTRODUCTION 


Recent investigations on the behavior of elastic slender rotor blades 
■undergoing relatively large deformations during operation, show that non- 
linear phenomena have considerable influence on this behavior. These non- 
linear phenomena are due to the inclusion of moderately large deformations 
in the elastic, inertial and aerodynamic operators associated with this 
problem. A detailed review of recent research on rotary ving aeroelasticity 
with on emphasis on the importance of moderately large deformations has been 
presented in Ref. 1 and is beyond the scope of this report. 

Recent emphasis on wind energy conversion, using large horizontal axiB 
wind turbines, employing two-bladed, hihgelesa rotor configurations having 
rotor diameters varying between 120 to 500 ft. have added a new stimulus 
to the study of large flexible, highly pretwisted blades. In a recent study 
by Friedmann (Ref. 2 ) it has been noted that efficient construction and opera- 
tion of wind turbines requires that the vibratory loads and stresses on the 
rotor itself and the combined rotor tower system be reduced to the lowest 
possible levels. Thus, structural dynamic end aeroelastic considerations 
are of primary importance for both the design of wind turbines and the compar- 
ison of various potential wind turbine configurations. 

In view of this need, it was felt that a careful, fundamental deriva- 
tion of the equations of motion for slender rotor blades, possibly highly 
pretwisted, undergoing relatively large deformations during their operation, 
is required. These equations could be used as a basis for future studies, 
into which nonisotropic material behavior, such as required for the treatment 
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of composite blades, could be easily incorporated. The main objective of 
the present study is the derivation of such a set of equations. 

A fundamental work in this field was that of Houbolt and Brooks (Ref. 

3) where equations of equilibrium for the coupled bending and torsion of 
twisted nonuniform blades were derived. Although some nonlinear effects were 
included in their derivation, their final results can be considered as a 
linear representation of the problem. Following this work, other researchers 
presented derivations of equations which include additional nonlinear terms. 
These include, for example, the work of Arcidiacono (Ref. 4), Friedmann and 
Tong (Ref. 5), Hodges, et al (Ref. 6-8) and a recent work by Friedmann (Ref. 
9). The most detailed and comprehensive derivation of a set Of nonlinear 
elastic equilibrium equations is presented by Hodges and Dowell (Ref. 8). 

There the equations are obtained by two complementary methods, Hamilton's 
principle and the Newtonian method. 

In the present work a set of nonlinear elastic equilibrium equations 
of a blade are presented. These equations are derived for the case of small 
strain and finite rotations. The derivation includes several assumptions 
which are presented during the presentation. The equations are developed with 
respect to two different systems of coordinates. In each case the derivation 
is done using two complementary methods ; the Newtonian method and the principle 
of virtual work. Finally, the equations obtained in the present study are 
compared with those obtained in the previous studies. 

This report is a modified and abbreviated version of Reference 13, 
which contains a considerable amount of additional details. Furthermore, 
it should be noted that this report represents the first in a series of 
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three reports which document the research which has been performed under 
the grant. The second report (Ref. 20) deals with the aeroelastic 
stability and response problem of an isolated horizontal axis wind turbine 
blade. This report also contains typical single-blade aeroelastic stability 
boundaries together with blade response studies at operating conditions for 
the MQD-0 wind turbine, currently in operation at NASA Lewis Research Center. 
The third report (Ref. 21) deals with the aeroelastic response and stability 
of a coupled rotor-tower configuration corresponding to the NASA/foOE 
Mod-0 machine. 
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2. BASIC ASSUMPTIONS 


The geometry of the problem is shown in Figures 1 through 3« The 
following assumptions will be used in deriving the equations of motion. 

1) The blade is cantilevered at the hub, the feathering axis of the 
blade is preconed by an angle 0^. 

2 ) The blade can bend in two mutually perpendicular directions 
normal to the elastic axis of the blade, and can also twist around 
the elastic axis. The boundary conditions are those of a canti- 
levered beam. 

3) The blade has an arbitrary amount of pretwist which is assumed to 
be built in about the elastic axis of the blade. 

4) The blade cross section is symmetrical about the major principal 
axis. It has four distinct points: 

i) Elastic Center (E.C.) - the intersection point between the 
Elastic Axis (E.A.) and the cross section of the blade 

II) Center of Mass (C.G. ) 

III) Tension Center (T.C.)-the intersection point between the 
Tension Axis (T.A. ) and the cross section of the blade 

IV ) The Aerodynamic Center (A.C.) 

As shown in Figure 3 the C.G - E.C offset is denoted by X^., the 
T.C - E.C offset is denoted by and the A.C - E.C offset is denoted 

by X A , where it is understood that the offsets shown in Figure 3 are 
considered to be positive. 

5 ) The strains in the blade are always small, but the rotations can 
be finite (for additional details see Appendix b). 
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3. ELASTIC EQUILIBRIUM EQUATIONS OF THE BLADE 


In this section, the equilibrium equations of the deformed blade are 
given. It is assumed that the blade can be considered to be a deformable, 
slender rod, made of linearly isotropic, homogeneous material. As formerly 
indicated, the analysis is restricted to the case of small strains and 
finite rotations. Appendix A gives a brief summary of sane well known rela- 
tions of nonlinear deformat ions. In Appendix B expressions for rotating and 
strains of a deformed slender rod are derived and the force and moment 
resultants are obtained. In Appendix C the equilibrium equations are derived 
systematically with respect to the deformed as well as the undeformed system 
of coordinates, using the Newtonian method. In Appendix D the same equations 
are derived using the principle of virtual work. 

In this study the Bemoulli-Euler hypothesis is assumed to apply. This 
hypothesis is usually stated as: "Plane cross sections which are normal to the 

elastic axis before deformation remain plane after deformation (except for 
negligible errors due to warping) and normal to the deformed axis." Furthermore, 
it is also assumed that strains within the cross section can be neglected, 
and the warping is very small so that its influence is negligible, besides 
its effect on the torsional stiffness. (For a more accurate approach other 
warping effects can be included as shown in Appendix B. ) 

As shown in Figures 2 and 3, before the deformation of the elastic 
axis of the blade, which is the line that connects the shear centers of the 
blade cross sections, coincides with the axis. The y Q axis is 
orthogonal to Xq and lies in a plane parallel to the hub plane, while z Q 
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is perpendicular to Xq and y Q . It is clear that Xq,Yq,Zq is a rectangu- 


lar Cartesian system. As shown in Figures 2 and 3, e ,e ,e are unit vectors 
- x y z 

in the directions x^,y 0 ,z Q , respectively. According to the Bemoulli-Euler 

hypothesis and the other accompanying assumptions, during the deformation the 

triad e^, e y , e^ is carried in a rigid form, composed of translation and - 

rotation, to the new orthogonal triad e',e',e' (shown in Figures 2 and 3). 

x y z 

The unit vector e ' is tangent to the deformed elastic axis, while e and 

x ^ y 

e are rotated around it to the position of e ' and e ' . 
z y z 

It is assumed that the blade is acted upon by a distributed load, p, 

per unit length of its undeformed axis, given in component form by: 


P = P e 
x x 


p e' + 

y y 


A, 

p e' 
*z z 


( 1 ) 


This load, p, includes body forces, surface tractions and inertial loading. 
There is also a distributed moment, q, per unit length of the undeformed 
axis of the rod, given by: 


* = Vi + Vy + 



( 2 ) 


This also includes body couples, moments of surface tractions and moments of 
inertial loading. The loads and couples, p and q, and their derivatives 
are assumed to be continuous. 

Then the exact equilibrium equations are obtained in Appendix C as 
Equations (C-7) and (C-8): 
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T + k (m + tM + q_) 

,x y v z,x y ^z 

- k (m - tM + q) + p = 0 
2 y>x z ^y *x 

-(M + k M + tM + q ) 

z,x z x y ,x 

+ kT-t(m + k m - tM + q) + p =0 
y y y x z t y 

(M + k M - TM + a ) 
y,x y x z y ,x 

* V - T(M z,x + *z\ + Tl ir + + J>z - 0 

M -km - KM+a«0 
x,x y y z z ^x 

The first three equations are basically from force equilibrium rela- 
tions in the e ' ,e ' , and e' directions, respectively. The fourth 
x y * 

equation is the moment equilibrium relation in the e^ direction. The equi- 
librium of moments in the directions e ' and e ' are also satisfied. 

y 2 

T is the axial tension in the blade. M ,M , and M are the compo- 

x y z * 

nents of the elastic moments, M is the torque, while M and M are the 

x y z 

bending moments, k and k are the curvatures, while t is the twist of 

y z 

the deformed elastic axis. The expressions for the components of the moments 
are obtained from Equation (B-43) of Appendix B: 

M = GJt 
x 

M y " " E3 23 K y * EI 33 IC z + tt oc 
M z " E3 22 K y + E3 2 5 \ “ ^oc 
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where E is young's modulus, G is the shear modulus, J is the torsional 
stiffness of the blade, and 811(1 £>j are the flexural moments of 

inertia of the cross section (effective in carrying tensile stresses ) around 


axes parallel to the directions e^ and e z which pass through the point 
(y QC ). This point, whose coordinates are (y^z^), * s the tensile 
center. The moments of inertia are given by Equation (B-U2 ), which is: 


*22 - JJ 

- V 2 % ^0 

A 


hi ' II 

( *oc - V <Z oc - *0 ,dy 0 ^0 

A 


hi * If 

(z oc * z 0 )2 ^0 4z 0 • 


A 


(5) 


Furthermore, it is assumed that the blade cross section is symmetric 
about the r\ axis (see Figure 3). The moments of inertia about r\ and an 
axis perpendicular to T), which pass through the tensile center, are 
denoted by 1^ and X> , respectively. Then (see Figure 3): 


hz * h cos ' "a * : 3 311,2 ®g J 

hi ' ( h - I 5 ) “ in ®G COS «G (6) 

hi ~ h 8ln2 ®g + h cos2 ®g • 

From Figure 3, the following relations follow: 


oc 


X II cos ®G 


OC 


X II Sin e G 


(7) 
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It is assumed that X is small enough so that the expressions 

T z and !IV in Equation (4) are at most of the magnitude of the 

oc oc 

other terms in the equation. 

Furthermore, it should be emphasized that Equations (4) were 

2 

obtained after neglecting terms of order £ compared to unity (for more 
details about the ordering scheme, see Appendices E and c). According to 
the assumptions of Appendix B, rotations are of order £. 

Within the order of approximation implied by neglecting terms of 

2 

order £ , compared to unity, the equations of equilibrium are simplified 
and are given in their final form by Equations (c-10): 


T +KM -KM + t(k M + k M ) 
,x y z,x z y,x y y z z 


K q 

y z 


k q + 
27 


M 

z,xx 


(k + tk )m 
z,x y X 


(t + k k )M - 2T M 

y Z y y,x 


+ k. T + K. q 

y z^x 


- q + p 
y 


M + (k - Tk )m - (t - k k )m - 2t M 
y,xx y,x z x ,x y z z z,x 


+ KT-Kq+q - Tq_ + P 

z y^c t,x ^z z 


M -KM -KM + q = 0 

x,x y y z z tc 
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In order to complete the formulation of the problem, the boundary- 
conditions must be also taken into account. For the present case of a canti- 
levered blade, the boundary condition at x^ = 0 corresponds to a clamped 
root and at the tip, = i (where i is the length of the blade), free 
end conditions apply. Thus 


for Xq - 0: 
for Xq = i: 


v = w = v 



w * ♦ = 0 
>x 


M = M = M 
x y z 


= 0 


( 9 ) 


V and V are the resultant shearing forces at the blade cross section 

y z 

and they are given by (see Equation (C-6) of Appendix c): 


V =-(m +KM+TM+q) 
y z,x z x y ts 


( 10 ) 


V 

z 


- M 


y>x 


+ K M 
Z X 




The underlined terms in Equation (lO) disappear in the case of a free edge 

when M = M - M =0, also, 
x y z 

The deflection of a point on the elastic axis is given by W where 
(see Equation (B-6) of Appendix B): 


W = ue + ve + we • 
x y z 


( 11 ) 
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The new triad, e'.e'.e', which is tangent to the deformed coordi- 
x y z 

nates of the blade is given by Equation (C-ll): 


e + S._e + S_,e 
x 12 y 13 z 

S *e + e — + S 0 ,e 
21 x y 23 z 

S,, e + S,„e + e 

31 x 32 y z 


where S. . are functions of w , v and <t> . When finite rotations are 
ij ,x ;X 

considered, these depend on the sequence of rotations which transform 

(e , e ,e ) to (e' ,e' r e' ), as one can see from Equations (B-10), (B-13) and 
x y z x y z 

(B-lM of Appendix B. If the sequence chosen is a finite rotation about e , 

z 

followed by rotation about e , followed by rotation about e , then, 

y x 

2 

neglecting terms of order £ compared to unity, the transformation given in 


Equation (12 ) is defined by Equation (B-13): 


e + v e + w e 
x ,x y ,x z 

- (v +<>w )e + e + t>e 

>x ,x x y z 

- (w -< > v)e-( <> +vw)e+e 

,x ,x x ,x ,x y z 


With relations (13) the curvatures and twist are then given by (B-16): 


v + w 
,xx ,xx 
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K = W - ©V 
Z ,XX , XX 

T = i + v w 

,x ,xx ,x 


(l*rt>) 

(l4c) 


The term v w in Eqs. (l3) and (l4) was introduced by Wempner 
» x ,x 

(Ref. 11 ) and was later shown to be significant for rotor blades by Kaza 
and Kvaternik (Refs. 15 and IT). Substitution of Equations (lk), ( 6 ) and 
(7) into Equation (b) implies: 


M y - -E(L, - I,) sin e G co.# G (v n t *»„) 

- E(I 2 e G + Ij cos 2 S G )(» iXjt - (v m ) + TX n sin e G 

M z = Edj cos 2 9 G + Ij sin 2 9 g )(v )TV + i « n ) 

+ e(i 2 - I 5 > sl ” ®G cos 6 a (w ,xx - - IX II cos 9 g 


• (15) 


Substitution of Equations (lk) and (l5) into Equation ( 8 ), using 

2 

Equations ( 6 ) and (7) and neglecting terms of order e compared to 
unity, implies : 


T ,x + \xx tEJ 22 <T ,xx + 'Kxx ) + E3 2 3 (v ,xx - 'Kxx 9 ,: 
+ w ,xx^ [EI 22 (v ,xx + 'Kxx 1 + EI 23 \xx 15 ,x 
+ ",xx lEI 23 (v ,xx + *“,xx ) + EI 33 ( ",xx * ■Kxx”, 

- t ,xx^ [EI 23 v ,xx + EI 33 (w ,xx ' 'K** )n ,x 

+ (v + 4w )q - (w - <^v )q_ + p =0 

,xx ,xx z ,xx ,xx 7 x 


( 16 a) 
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- [EL 0 (v + ) + EL,(w - ) - Ty ] 

z2 ,xx ,xx ~3 ,xx ,xx J oc ,xx 

- GJ(* w +vww - *<> v ) + EL, * v 

,x ,xxx ,xx ,x ,xxx ,x ,xxx z3 >XX , XX 

2 

+ EI,,(* w + v ww - v +2vw ) 

33 >xx ,XX ,XXX ,X ,XX ,XX ,XX ,XX , XX 

+ 2(* + v v )[EL, v + EI,,(w - <V ) - Tz ] 

,x ,xx ,x ^3 ,xx 33 v ,xx , xx oc 

+ (v + - z * )T + (w - ♦v )q 

— — — ,XX ,XX TC 


,xx 


,xx oc ,xx 


- (♦ +v w)q-q_ +p 
>x >xx ,x V ^z,x *3 


(l 6 b) 


xx 


- [EL,(v + ) + EI,,(w - *v ) - Tz 1 

L3 , xx ,xx 33 >xx , xx oc , 

+ GJ( (|, v + ^ w + v v w ) - EL, * w 

,x ,xxx ,x ,xxx ,xx ,xxx ,x ^3 ,xx ,xx 

- EL„(* v + ‘* 4> w + v v w) 

c.2 ,xx ,xx ,xx ,xx ,xx ,xxx ,x 

- 2( *,x + v ,«»,x )lE W T , 1 ot + *",** ) * EI 23 w ,«‘ ^oc 1 ,: 

+ (w - "V + y * )T - (v + V )q 

,XX ,XX OC ,XX ,XX ,XX 

- (♦ x + ▼ _y _)q_ + Sr x + p * ■ 0 

; X f XX jX Z J ;X 2 


(l6c) 


[Gj(t> + v v )] + 

,X , XX ,x ,x 


■ w + 4$ v w ) 

r XX ,XX ,XX , XX 


EI 23 (v ,: 

+ < EI 33 - EI 22 )t ’ r ,xx w ,xx + 

* Tly o<>=oc - - z oc (v ,x* + *",>=x )1 + \ ■ 0 


- V )] 
, XX , XX 


(l 6 d) 


Hie system of Equations (l 6 ) contains four equations with four 

unknowns being represented by v, w, ♦ and T. The boundary conditions 

(9) together with Equations (10) and (15) become: 

for x~ ■ 0 : v = w = v =* w = ♦ = 0 , (l7a) 

u »x >x 

for Xq = i: T ■ 0 , (l7b) 

and 
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-V _ \ = [E(L cos 2 9 + I sin 2 ©~)(v + <5> w ) 

y(M =M =0) z G 3 G ,xx , xx 

x y 

+ Edj - Ij) sin 6 G cos e G („ >xx - <V )Xx )-T X n cos e G ] )X 

+ q = 0 (17c) 

z 

V m -r.'i = -(E(L - I,) sin 9 cos ©„(v + V ) 

z(M=M=0) z 3 G G ,xx , xx 

X z 

+ E(x, Sin 2 e + I cos 2 0 )( w - ‘•V )-T X sin © ] 

£ G 3 G ,xx ,xx II G ,x 

+ = 0 (I7d) 

M = G J + v w) = 0 (I7e) 

X ,x ,xx ,x 

-M y - Ktt, - I 3 ) sin e G cos S p (V (3Dt+ *„) 

+ E(lg sin 2 9, + Ij cos 2 e Q )(W (XX - V xx ) - 0 (I7f ) 

M z - cos 2 e Q + Ij sin 2 »„)(▼„ + *» )XX > 

+ Edj - Ij) sin e G cos e G (w jia[ - «v >xx ) = 0 . (17g) 


The equations of equilibrium presented above were derived in the 

directions of the deformed coordinates. As pointed out in Appendix C, the 

equations can also be derived in the directions of the undeformed coordinates 

e ,e ,e . In this case the distributed force, given previously by Equation 
x y 2 

a), can be taken in the form: 


P 


p e + 
x x 


p e + 

y y 


(18) 
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while from similar considerations the distributed moment can be written as: 


<1 


Vr + + 


t y 


q e 

2 Z 


(19) 


The equilibrium equations are given by (C-24 )„ In the present case 

the triad e ' ,e ' >e ' is given by Equation (13). Then Equations (12 ) and 
x y z 

(13 ) imply: 


S „ = v 
12 >x 


S 21 = - (v ,x + 






- ); 



( 20 ) 


- (<t> + 



Substitution of Equation (20 ) into Equation ( C-2 ^ ) and neglecting 
2 

terms of order £ compared to unity, yields the following equilibrium 
equations : 


(T + (v + )M - (w - ♦v )M + w * M 
,x , x z,x ,x , x y,x ,x ,x z 

+ v <tl M+(vw -v w)m] 

,x ,x y ,x ,xx ,xx ,x X ,x 





0 


(21a) 


-[M +w M + (♦ +v w)M+( l)> +vw)M ] 

z,x ,xx X ,x ,xx , x y ,x ,x y,x ,x 


+ 'V’.rV* + 


0 


(21b) 


[M +v M-( 0 + v v )M - ♦M ] 

y,x ,xx x ,x ,xx ,x z z,x ,x 

+ (w ,x T) ,x- ^.x^x’.x* Vx + *z ’ 0 


(21c) 
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M - (v + ‘V )M - (w - <*v )m 

x,x , xx , xx y , xx , xx z 


+ V 


o+w q+q = 0 

,x t ,x x 


(2 Id) 


Substitution of Equations (15 ) into Equations (2l), using Equations 

2 

(6) and (7), and neglecting terms of order e when compared to unity, 
yields the following equations 


(T + v [EL„(v + *w )+ Ei,,(w - *v )] 

,x , xx , xx , xx , xx ,x 


+ ) + El,, (w - <V )] 

xx , xx 33 , xx ,xx ,x 


+ w ,x [EI 23 (v , 

+ 1 ',x'* [EI 22 (v ,xx + " u ,xx ) ' EI 2 3 “,xx 1! ,x 

- v ,xt* [EI 23 \xx + EI 33 (V ,XX - \xx>”,x 

+ GJ(0 + v w )(v w - v w )) 

,x ,XX , X , X ,XX ,X3C ,x ,x 

+ (v q) - (w q) +p = 0 

,x z ,x ,x T >X 


(22a) 


- [E V T ,xx + ^,xx J + EI 23 ( ”,xx - 2 * v ,xx ,] . 


XX 


- [GJ(0 + V w )w ] 

,X ,xx ,x ,xx ,x 

+ [EI-.,(0w + v w w - £ v )] 

33 ,xx ,x ,x ,xx ,xx , xx 

- (El,, V W 2 ) + [v T + (Ty ) - (*Tz ) ] 

33 ,x ,xx ,x ,x J oc ,x oc ,x ,x 


- S,x + <W ,x + ?y - 0 


(22b) 


.xx 


- [EI_,,(v + 20 w ) + El,, (w " - 0v )] 

^3 , xx ,xx 33 ,xx ,xx , 

+ [GJ(0 + V w )v 3 - [ELo 0(v + )] 

,X ,XX ,X ,XX ,X 1 ^2 v ,xx ,xx ,x 

- (EL_ V 2 w ) + [w T + (Tz ) + (OTy ) ] 

22 ,xx ,x ,x ,x oc ,x oc ,x ,x 
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(V ,X Sc^X + VjX 


(22c) 


+ p = 0 

*z 


[GJ( •** + v w )] + EX,^(v 2 - w 2 + M> v w ) 

,X , XX ,x ,x , xx , XX , XX ,, — 


i XX 


OC , XX 



)[v w 

+ 4 >(w 2 

2 

- V 

:2 

, XX , XX 

,xx 

y 

+ 

1 

> 

I 

* 

•V 


,xx 

OC , XX 

y 


XX 


XX 


)] 


+ v q + w q+q 
»x 7 ,x ’z at 


= 0 


(22 d) 


The boundary conditions remain the same as in Equation (17) and it 
is only required to write q z and as functions of q^,q^, and <l z * 
According to Equations (C-l6) and (C-ll): 


Sr - V S 2lSc + S 23\ ; 

_ - . (23) 

S = S + S 31 Sr + s 32 Sr • 


The equilibrium equations in Appendix C were derived by the Newtonian 
method. In Appendix D the two sets of equations, with respect to the 
two different systems of coordinates, are derived using the principle of 
virtual work. The equilibrium equations which are obtained from this 
procedure are identical to those obtained in Appendix C, within the 
approximations inherent in the present theory. One of the advantages 
of the second method is that it also provides the appropriate set of 
boundary conditions, which is sometimes difficult to obtain using the 
Newtonian method. It is shown that the boundary conditions of the blade, 
as stated in Equations (9) or (17), are in agreement with the boundary 
conditions obtained by the second method. 
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The equations of equilibrium can be further simplified by taking 
into account some common properties pertaining to helicopter and vind 
turbine blades. These blades are usually stiffer in lag than in the 
flap-wise direction, thus: 

EL, > EIj , GJ . (24a) 

T _ 

The geometric pitch angle 0_ has an absolute value less than 

Or 

45°. Therefore, according to Equation (6) 

EIgg > EI 33 ,GJ . (24b) 

Using Equation (24b) together with the ordering scheme (e neglected 
compared to unity) enables one to neglect a considerable number of addi- 
tional small terms; the resulting equations are given below: 


< T + V [EI 22 ( LS + 


XX 




XX ,x 


- *v )] 
,XX .XX ,x 


+ ".x^V^g + !L S£ ) + EI 53t_ 

4 “.x ;t(E ^ (Y ,xx* ^'■xx ) + E1 23 v .xx 1! .x 

- V (^[EL, v + El,* w ]} + GJ 41 (v w - v w )) 

.x 1 - ^3 .xx 35 ,xx »x ,x ,x .xx , xx ,x ,x 

+ ^x^U- ( “,x V,x + k 


(25a) 


[EI 22 (t ,xx* \xx } + E V w ,xx - - EI 33 \xx ] ,xx 

+ [v ,x T+ ^oc>,x- ( ***oe>,x 1 ,x 


- Vx* ( “,x ^,x + Py = 0 


(25b) 
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, XX 


- [EL,(v + 20 v )+ EI„(w - 1 V )] 

z3 , xx , xx 33 ,xx ,xx , 

+ [GJ^ + v w )v ] - [EL, <Kv + "Hr' )] 

- (EL, V 2 w ) + [w T + (Tz ) + ( *Ty ) ] 

<22 /XX ,x ,X ,X OC ', X J OC ,X /X 

. ~ N, ~ ~ 

- <v ,*«* ) ,* + V,* + p * = 0 


XX 


(25c) 


[GJ($ + v w )] + EL,(v 2 - w 2 + 4$ v w ) 

,x ,xx ,x ,x ^3 ,xx ,XX , XX , XX 


< I Sj- E %* )[T ,J , ,a + * (w !xx- v ! 


XX 


.)] 


- T[z (v + 'V ) - y (w - “hr )] 

OC , XX , XX J OC v , XX /XX 


+ V 


q_ + w q_ + q_ 
,x ,x jc 


= O 


( 25 d) 


Usually, in the case of rotating blades large tensile forces are caused 
by the centrifugal forces; therefore, the nonlinear contributions of the 
bending and torsional moments to the tensile force are very small. In 
this case it seems to be justified to keep only the principal terms of 
this contribution (twice underlined in Equation (25a)) and neglect all 
the other terms associated with it (once underlined in Equation ( 25 a)). 
In fact, it seems that neglecting the twice underlined terms will also 
not affect the results in a significant manner. 


A further simplification can be obtained if all stiffhesses are 
approximately of the same order of magnitude, which means : 


*22 

I 33 



22 s 


(0.5 - 2) 


In this case. Equations (22 ) turn out to be : 


(26) 
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ft* ^.xx 1 * E W».xx - 

* ^ * Br 53 ( L« ' * t .xx ) \x 

* ».xftft^ \xx* 5a W xx )] .x - ZJH^I y -xx* EI 33 »„>!., 

+ aJ *.* (T .«*.« 1 r ,xx v .x^.x * <v 5 x ) ,x - (v ,x 5 y ) ,x + *x'° 


- ^\xx * * K ,xx ) * E3 23 (v ,xx - 2 * T ,xx ) - EI ?? \xx J , 

- (0J *,x ' , ,xx ) ,x + [ \x T + ^oc\x - ( * T V ) ,x ! ,x 

- Vx * (k ,x \\x * p y ■ 0 

- XX + 2 ’",XX } * EI J3 ( *,xx - ‘ T ,xx ) * E1 Z2 *W, 
* <**,« v ,xx>,x * ft* T+ (ft oc\x * ^oc^.x 

- (t ,x \\x * V.x + *s " 0 


XX 


XX 


lGJ ft,x * T ,xx v ,x ,1 ,x * EI 23 (v !xx - ^ + ft 1 ?? - EI 22 ^.xx »,xx 
- Tl V (r ,xx * * l, ,xx ) - y <x>xx - ^.xx* 5 
* y ,x V V \ * *x * 0 


Hie underlined terms in Squat ion (27a) have the same meaning as 
those in Equation (25a). 

To facilitate the use of these equations for rotor-dynamics 
applications and to also simplify comparison of the equations in this 
revised version of the report with the previous version, the Equations 
(27a-d) are rewritten below using the principal moments of inertia of 
the cross section (Equations (6) and (7))» It should be noted that in 


(27a) 


(27b) 


(27c) 


(27d) 
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the previous version of the report, few terns in Equations (22b - d) were 
missing due to an algebraic error. In most copies of the 

report these terms were added in handwriting. 


T ,x + - 0 


(27aa) 


cos 2 0 + I, Sin 2 © )(v + Ow ) 

Cj ^ m y XX y XX 


.XX 


- Mh 

+ E(lg - Ij) sin 9 g cos 6 j(w - 2 * v 

- E < 1 2 ^ e G * V^xx^xx 

* ( v ,x T + [m ii (cos e G- * (, ,A ) , I , 'r 


(GJO w ) 

,X ,XX ,X 


- 0 (27bb) 


- [E(I 2 ' V Sl " ®G COS 9 G ( \xx + Zt “.xx* 
* E(3 2 3l “ 2 V T ) cos2 e G )(u ,xx - *' , ,xx ) 
+ E <% 


cos e + I, sin © )<*v _] + (GJO v ) 

G y vT . y XX y XX y X y XX y X 


* (\x T+ [' E! Il( si ” e G + * COS e G )1 ,x ! ,x 


- (v „ q) 


,x Sc^x + 


+ p„ 


(27cc) 


[GJ( 0 + V w )] +(!_-: I ) sin 0 cos 6 (v 2 - w 2 ) 

,x ,xx ,x ,x 72 3 G G ,xx ,xx 


+ (ij - l 2 )(cos‘ : 


9 - sin 0_ )v w 

G G ,xx ,xx 


- TX TT [sin 0_(v + ) - cos ©_(w - V )] 

II G ,xx ,xx G ,xx ,xx 


+ v 


x q v +W x q z +q x 

y JL j > A <u X 


= 0 


( 27 dd) 
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4. COMPARISON BETWEEN THE PRESENT ELASTIC EQUATIONS 
AND THE EQUATIONS OF OTHER STUDIES 


Elastic equilibrium equations for a rotor blade were derived by 
different researchers during the past twenty years, as was shown in the 
introduction. In this chapter a comparison will be made between the 
present derivation and some of the previous ones. For the sake of 
brevity these comparisons are concise, much more detailed comparisons 
can be found in Reference 18. 

4.1 Comparison with Houbolt and Brooks (Ref. 3) 

A set of equations equivalent to those of Reference 3 can be 
obtained from Equations (21a - d) by neglecting the nonlinear terms assoc- 
iated with the elastic moments. Performing these operations and replac- 
ing the moments by the appropriate expressions, as shown in Equation (15) 
results in the equations given below. It should be noted that nonlinear 
terms containing the displacements have been neglected in these equations 


P + (ql v ) - (q, w ) + p - 0 

,x ,x ,x y >x ,x *x 


, (28a) 


tE(3 2 COs2 ®G + b Sln2 ®G )v ,xx 

+ £(^-13) sin e G cos e G w #xx -TX II cos © G ] >xx (28b) 

^^x^x + q z,x " (w ,x ^x\x P y ° ’ 

- [ E(^ " I 5 ) sin © G cos e G V ^ 

2 2 

+ E(X, sin e,, + I, cos Q„}w -TX tt sin ©„] 

' £ G 3 G ,xx II G ,xx 


+ (Tw )+ q - (v q ) + p - 0 

,x y,x ,x x,x z 


(28c) 
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[GJ 4 J 


+ 

>x 


q_ w + 
^ ,x 


rKi 




(28d) 


If Equations (28) in this report are compared with Equations (15) 
and (l8-20) of Reference 3, the following observations can be made: 

1) Equation (28a) is identical to Equation (15) of Houbolt and 
Brooks, except for the terms (? v ) and (q" w ) 

« yX yX J ,X y X 

which are not present in Equation (15). 

2 ) Comparing Equation (28b) of the present study with Equation (20 ) 
of Reference 3, it follows that in addition to the terms con- 
tained in Equation (28b), Houbolt and Brooks' equation contains 
an additional term involving P , probably resulting from the 
assumption Q' 1 1 ® a . 

Furthermore, the term [ TX^ ^ sin © Q ] ^ which appears 
in Reference 3 is a physically nonlinear term. This term, which 
appears in Equation (22b) of the present study is associated 
with the term (<J> M^) ^ in Equation (21b) of the present study. 
If this term is retained, all other terms of the same order 
should also be retained. However, this was not done in the 
equations presented by Houbolt and Brooks. 

In the loading terms, the term, - (w q) , which 

)X A jX 

appears in Equation (28b) of the present study does not appear 
in Equation (20 ) of Reference 3. 

3) The comparison between Equation (28c) of the present study and 
Equation (19) of Houbolt and Brooks is analogous to the compar- 
ison given in the previous section, and will not be repeated 
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here 


4) Comparison of Equation (28d) and Equation (l8) of Reference 
3, shows that except for the terms which appear in Equa- 
tion (28d), additional terms containing P appear in 
those of Houbolt and Brooks' equation. These, again, are 
related to the assumption cr^ » a . Houbolt and Brooks 

also retain the term Tk i (k being the radius of 

& 8 . 

gyration of the cross sectional area). However, this term 

should be neglected within the assumption that strain is 

negligible compared to unity. Similar to what has been 

pointed out already in Item (2 ) of this comparison, from the 

nonlinear terms - v M.. and - w M in Equation 

, xx y , xx z 

(2 Id) of the present study, Houbolt and Brooks retain only 
the terms - T X TT sin ©_ v and TX_ t cos ©_ w , 

Or ^ XX 11 G f 3CX 

while apparently neglecting other terms of the same order. 

5) As pointed out in (2 ) and (4) above, it appears that 

Houbolt and Brooks assumed that in the expressions for the 

moments M z and M y (Eq. (4) of the present study), the 

terms Tz and fy cure much larger than the other 
oc oc 

terms. This seems to imply that the offset between the 
shear center and the tension center is the main contri- 
butor to the bending moments in the blade . This is a very 
special case which may be of limited importance from an 
engineering point of view. The assumption of the present 
study, that these terms are,, at most, of the same 
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magnitude as that of the other terms, seems to be more 
realistic. 

k.2 Comparison with Hodges and Dowell (Ref. 8) 

1) Comparison of the expressions for shows that the 

strain at the elastic axis, and the contributions due to 

bending all over the cross section, are exactly the same. 

The expressions for warping are different. Hodges and 

Dowell also add the terms (t) 2 + )©„ <J> and 

G,x ,x 

((t) 2 + 2 • The first expression is also present in 

the derivations of Houbolt and Brooks. The second term, as 

was pointed out by Hodges and Dowell thermselves, is neglect- 

2 

ed within the approximation that terms of order £ are 
negligible compared to unity. They retained this term for 
the case of very large torsional deformations, which imply 
that y <l> and | are of order £ . The occurrence 

of such large elastic twist is unusual for most wind turbine 
or helicopter blades and is not treated in the present work. 

It should also be mentioned that the case when squares of 
strains are not negligible, in comparison to the strains 
themselves, is very special, and it appears that for this case 
a more refined theory than the one presented in this study 
will be required. 

2) Comparing the shearing strains e _^ and as derived 

from Equation (B-25) of the present study with Equations 
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(25, 26) of Hodges and Dowell, shows some discrepancy which 
appears to be related to the previously discussed differ- 
ences in the warping function, and the expression for the 
twist. Furthermore, terms which represent products of the 
warping and curvatures in this work, are neglected by Hodges 
and Dowell in theirs. This neglect seems to be justified 
for the case of slender blades having a closed cross section. 

After obtaining the strain components, Hodges and Dowell derived 
the equations of equilibrium by two complementary methods, the Hamil- 
tonian principle (similar to what is done in Appendix D of this report) 
and the Newtonian method (similar to what is done in Appendix C of this 
report). It is obvious that different strain expressions, and different 
transformation relations between (e^ , e y , e^) and (e^ , e^ , e^), will 
yield different final equations. 
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5. CONCLUDING REMARKS 


A system of consistent, nonlinear, equations of equilibrium of a 

pretwisted wind turbine or helicopter blade which undergoes moderate 

deformation was systematically derived. The derivation contains, in 

addition to the basic assumptions listed in Section 2, some additional 

assumptions which are gradually introduced in the course of derivation. 

For the sake of completeness, these additional assumptions are briefly 

summarized. The blade is slender and its undeformed elastic axis is 

straight, the blade is made of elastic isotropic material. The Euler- 

Bernoulli assumptions are valid (for details see Appendix B) and warp- 

* 

ing of the cross sections due to torsion is neglected. Axial forces 
in the blade contribute to the bending moments, due to the offset be- 
tween the elastic center of the cross section and the tensile center. 

It is assumed that this offset is sufficiently small such that the mag- 
nitude of this contribution is, at most, of the magnitude of the other 
contributions to the bending moments (e.g., see Equations (B^3)). The 
strains in the blade are always small (less then 0.01 ), while the 

slopes due to elastic rotations are of order of magnitude e where 

2 

e « 0.2; furthermore, terms of order of magnitude e are neglected 
when compared to unity. Finally, it is assumed that deformations are 
changing gradually along the span of the blade, which implies that a 
modal expansion representing blade deformations would be restricted to 

•fr 

Except for the warping contribution to torsional rigidity. 
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the lower modes. 


Nonlinear structural problems require careful distinction between 
undeformed and deformed systems of coordinates for representing blade 
deformations. Therefore, in this study the final equations of equili- 
brium are presented in both the undeformed and deformed system. The 
general load components are also defined with respect to each of the 
systems . 

The orthogonal system of coordinates : x, y, z, used in deriving 
the equations of equilibrium in this study was found to be slightly more 
convenient than the curvilinear nonorthogonal x, q, £ coordinate 
system used in References 3, 6, 8, 9 and 17- The main advantage being 
a somewhat simpler derivation and slightly simpler final equations. 
Additional information on this topic is provided in Reference 18. 

An ordering scheme, such as used in this study, can simplify the 
equations considerably. The equations can be further simplified for 
certain blade geometries. It should be noted that the loading terms in 
the equations (forces and moments) were presented in a general form, 
without any approximations. Substitution of explicit expressions for the 
loading terms, and the application of an ordering scheme, enables one to 
identify and neglect a considerable number of additional small terms. 

Since their derivation, these equations have been used extensively 
in a variety of aeroelastic stability and response problems as indicated 
below: 

l) Calculation of coupled flap-lag-torsional aeroelastic stability 
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of hingeless rotor blades in forward flight (Ref. 19)* 

2 ) Aeroelastic stability and response calculations for an iso- 
lated horizontal axis wind turbine blade (Ref. 20 ). It should 
be noted that in this study, dynamic blade root bending 
moments were also calculated and found to be in satisfactory 
agreement with the loads measured on the NASA/DGE Mod-0 
machine. 

3) Aeroelastic stability and response calculation of a coupled 
rotor/tower horizontal axis wind turbine, simulating the 
behavior of the NASA/DCE Mod-0 machine (Ref. 2l). 

Finally, it is important to note that the equations derived in this 
study were used to investigate the large deformations of a cantilevered 
beam loaded by a concentrated transverse load at the free end (Ref. 22). 
The numerical results obtained were in very good agreement with experi- 
mental results , which indicates that these equations are reliable and can 
be used with confidence in a variety of applications. 
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Figure la. General Description of Helicopter Rotor Geometry. 
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Fig. 2b. Geometry of the Elastic Axis of the Deformed Blade 
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POSITION OF BLADE CROSS SECTION 
BEFORE DEFORMATION 


APPENDIX A 


DEFORMATIONS 

A detailed development of the expressions in this Appendix can be 
found in many books on elasticity (for example, Refs. 10, 11 ); therefore, 
these are only briefly repeated here for the sake of clarity. 

Consider a material point P in an elastic body, where the posi- 
tion before deformation is given by the position vector r (Figure A-l). 
The position vector r is a function of three coordinates, such that: 

r - r(x 0 ,y 0 ,z 0 ) . (A-l) 

The coordinate system shown in Figure A-l is an orthogonal Cartesian 
system with the unit vectors e x , e^, in the directions x D <yo ,Z 0 > 

respectively. Thus : 

7 Vi + Vy + Vi * < A ' 2 > 

After the deformation, the material particle is located at point 
P' (see Figure A-l) defined by the position vector R. If the initial 
coordinates of the particle are used as independent variables, then: 

R - R(x^ , y Q , z Q , t) , (A- 5) 

where R is also a function of time because the position of the point is 
a function of time (x^y^, z^ can be considered to be the coordinates of 


36 



the point at t = t Q ). 

If V denotes the displacement of the point, then: 


R » r + V . 


The base vectors of the point before deformation are defined as : 


g x =r ,x ; *y =r 


; g « r 
,y 6 z ,z 


After the deformation, the base vectors are: 


G = R ; G »R ; G “ R 
x ,x y ,y ’ z ,z 


The strain components are given by Equation (2-20 ) of Wempner in 
Reference 11 (where XQ,y 0 ,z Q is an orthogonal system): 


^ 55 | (G • G - l) ; 

xx 2- ' x x - - - - 


e = 
yy 


s i t* 


2 { \ ' % - X) » 


= 1 ^ 


zz 


? • 5, - 1) i 


: ■ e • i (5 • 5 ) J 

xy yx 2 x y 

i = e = ~ (5 • G ) ; 

xz zx 2 x z ’ 


s - e . - ~ (g • G ) . 

yz zy 2 y z 


(A-4) 


(A-5) 


(A-6) 


(A- 7) 
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Fig. A1. Position of a Material Particle Before and After the Deformation 
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APPENDIX B 


MECHANICS OF A DEFORMED ROD WHICH IS SLENDER AND 
STRAIGHT BEFORE DEFORMATION 


B.l General Expressions 

A straight slender rod is shown in Figure B-l. Every material point 
in this rod is described by a rectangular Cartesian system of coordinates, 
V y 0 ,z 0 . The coordinate x^ is identical with the elastic Mis of the 
rod, defined as the line which connects the shear centers of the cross 
sections of the rod. It is assumed that the elastic axis is a straight line. 
In this case, denotes length along the elastic axis of the undeformed 

rod, while y 0 and z^ denote lengths along lines orthogonal to the unde- 
formed elastic axis. 

Before the deformation the position vector of every material point is 
given by: 



while after the deformation, at time t, the new position vector is : 

5 53 *(*<) » y 0 > z 0 » • (B-2) 

The displacement of the particle is : 

V = R - r . (B-3) 
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Looking at a particle, which before deformation lies on the elastic 


axis, its initial position vector is: 



and its position after the deformation is given by: 


(B-4) 


q R - 5(xjj,0,0,t) . (B-5) 

The displacement of this point is denoted by: 

W = V(^ , 0 , 0) » ue x + ve y + we z . (B-6) 

The base vectors of the unde formed rod are simply (according to Eq. 
(A-5 )) the orthonormal triad: 


g ■ e 
x x 


5 ®y ~ e y 5 g z =e z ' 


(B-7) 


while the base vectors of the deformed rod are (according to (A-6)): 


G *5 =* (r + V) = e + V 

X ,X ,x x ,x 


G =* R ■ (r + v) = e + V 
y >y >y y >y 


(B-8) 


G 


R = (r+v) = e + V 


and at the elastic axis, a set E ,E ,E is defined as : 

’ x’ y* z 




The strains at any point are calculated by using Equations (A-7). 

The motion which carries the rectangular lines of the undeformed rod 

into the curved lines of the deformed rod, carries the initial tangent unit 

vectors e ,e ,e to the current tangent base vectors, G ,G ,G , 
x y z x y z 

respectively. This motion can be looked upon as two successive motions: 

First, the triad e ,e ,e is rigidly transformed and rotated to the 

x y z 

orientation of an intermediate orthonormal triad e',e',e'. Next, the inter- 

x y z 

mediate triad is deformed to the triad G ,G ,G which means changing the 

x 7 y z * 0*0 

angles between the vectors as well as the length of the vectors. The 
procedure which was described above is illustrated by Figure B-2. 

Consider a triad e ,e ,e which is positioned on the elastic axis 

x 7 y z •• - 

of the rod before deformation. In the stage of rigid transformation and 

rotation this triad is carried to the triad e' >e' ,e' , respectively. With- 

x y z 

out any loss of generality, let us assume that e^ is carried in this 
stage to the direction of E^, which means that it is tangent to the 
elastic axis of the rod after deformation. 

If the rotation of the triad e ,e ,e to the position e , ,e , ,e' 

x 7 y z x 7 y 7 z 


hi 



is relatively large, it cannot be described by a vector and it is treated by 


means of Euler angles (see, for example, Novozhilov, Ref. 12, Chapter VI, 

p. 205). If it is described as a finite rotation 0 about e , followed 

z z 

by a rotation 9 about e , followed by a rotation 0 about e , then 
J y y J x x’ 

the triad e' .e’,e' obtained after these three rotations, is given by: 
x y z 


(cos 9 cos 9 )e + (cos 9 sin 9 )e 
y z x y z y 


- sin 9 e 

y z 


, (B-lOa) 


e ' = (sin 9 sin 9 cos 9 - cos 9 sin 9 )e 

y x y z x z x 


+ (cos 9 cos 9 + sin 9 sin 9 sin 9 )e 
x z x y z' y 


+ (sin 9 cos 9 )e 
x y z 


, (B-lOb) 


(cos 9 sin 9 cos 9 + sin 9 sin 9 )e 
x y z x z x 


- (sin 9 cos 9 - cos 9 sin 9 sin 9 )e 
x z x y z y 


+ (cos 9 cos 9 )e 
x y z 


. (B-lOc) 


Equation (B-IO) are identical to Equations (A2 ) of Hodges and Dowell (Ref. 
8), after replacing 9 x , 9^ and 9 ^ by 9, and X, respectively. 

Consider the deformation of an element dx^ on the elastic axis 
of the rod, as described in Figure B-3. The procedure is as follows: 
First, the element is carried in a rigid body translation that does not 
appear in Figure B-3. Then the element is stretched by an amount 



6 about 

Z 


u x dx^ to the position A - (T) . Then the element is rotated by 
e z while point (l) moves a distance v x dx^ to location (2) , followed 
by a rotation -© , while the element tip moves from location (f) to @ • 

y 

Finally, the element in its position A - (j) is rotated by an amount © x 
about itself. From Figure B-3, the following relations are obtained: 


sin © 


1 


2 2 2 

l+2u +u +v +w 

,x ,x ,x ,x 


, (B-lla) 


cos © 




2 2 
1 + 2u + u + v 
t* i2L 


Jl + 2u + u 2 + v 2 + w 2 
V ,x ,x ,x ,: 


(B-llb ) 


8 in 0 « 

z 




2 2 

1 + 2u + u + V 

,x ,x ,x 


(B-llc) 


cos e 

Z 


1+ u 


■ 

V 


2 2 
1 + 2u + u + v 

,x ,x ,x 


(B-lld) 


The quantity u is of the magnitude of strain, as will be shown 

later. Assuming that the strains are small (say, £^ < 0.01 ) which is the 

case for most engineering materials, then the expression u will be 

neglected in comparison to unity. Next, it is assumed that v ,w and 

,x ,x 

are quantities of magnitude equal to or less than £ (in our case, 

2 

£ = 0.2 ), and quantities of the magnitude of £ are negligible compared 


to unity 



Using these assumptions and expanding sin © x and cos © x into 
series, one obtains: 


sin © 


e 


sin 9 = -w ; sin 9 = v ; 

y ,x z ,x ’ 


cos 


9 = 1 


cos 


e = 


1 ; cos 9=1 . 

z 


(B-12 ) 


At this stage, in order to be consistent with the usual notation in 
the literature (for example, Reference 9 ), ® x is replaced by Sub- 

stituting Equations (B-12 ) into Equations (B-10) implies: 


e + v e + w e 
x ,x y ,x z 


e ’ = - (v + <w )e + e + 

y ,x ,x x y z 


(B-13) 




-(w - )e - (0 + v w )e+e 

,x ,x x , x ,x y z 


Because the rotations are treated an finite, it is not surprising 


that the triad e*,e' ,e' depends on the sequence of rotations. If the 
x y z 

sequence consists of rotation © about e , followed by a rotation © 

y y ^ 

about e , and finally, a rotation © about e , as described in 

Z XX 

Figure B-4, then the expressions become: 


e' 


e + v e + w e 
x ,x y ,x z 


e ' ■ - (v + <>w )e + e + (♦ - v w )e 

y ,x ,x x y ,x ,x z 


-(w - <JV )e - ♦£ + e 

,x ,x x y z 


(B-14) 



The triad (B-l4) differs from (B-13) by terms of second order. Other 
different sequences of rotation will yield other triads which will differ 
from each other by second order terms. Therefore, it is most Important to 
retain one particular triad during a complete derivation, and consistency 
with this selected triad. In the theory of space curves there are three 
important quantities, defined as (for example, Wempner , Ref. 11, Eqs. 

(8-19) - (8-21)): 


K 

y 


e' 


x,x 



y>x 


K 

Z 




A 

e 


» 


z,x 


> * 


T 




(B-15) 


where and are curvatures, while t. is the twist. Equation (B-15 ) 

represents the exact expressions for the curvature and twist when strains 

are neglected compared to one. If e', e' , and e' are given by Equation 

x y z 

(b- 13), then: 


e + w 


e 


x,x 

,xx y , xx 

z 


e ' = 

-(v + $ w 

+ * w )e 

+ ♦ e 

y»x 

,xx ,x ,x 

,XX X 

,x z 

^ i 

e = 

- (w - 0 V 

ov )e - 

(♦ +v w +vw )e 

z,x 

,xx ,x ,x 

, XX x 

,x ,xx ,x ,x ,xx y 


(B-l6a) 


Substituting expressions (B-l6a) into Equations (B-15), and assuming 
that v , w , and are of the same magnitude, and neglecting 

y XX y XX y X 

2 

again terms of the magnitude of e compared to unity, implies: 
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(B-l6b) 


K = V + ♦w 

y ,xx , xx 

K = w - *V 

z ,xx ,xx 


T = ♦ + V W 

,x ,xx ,x 


The term v w in the expression for twist is a well known 

} XX t X 

term in the theory of rods (see for example, Reference 11, p. 390, 

Eq. 8-152d). It has also been used in rotor dynamics by 
Kaza and Kvatemik (Ref. 15). 


Fran the definitions (B-15 ) and the orthonormality conditions of 


x y ’ z J 


it is 

clear 

that: 

e ' 

= K 

A, A, 

e’ + k e ' 

x,x 

y 

y z z 

e ' 

= -K 

A A, 

e* + t e ’ 

y>x 

y 

x z 

e ’ 

= -K 

e' - t e ' 

z,x 

Z 

x y 


(B-17) 


which can be easily verified by substitution of expression (B-17) into 
the definitions (B-15). 


B.2 Bernoulli-Euler Hypothesis 

At this stage, it is necessary to find an expression for R. This 
always requires certain assumptions. In the present case, the well known 
Bernoulli-Euler hypothesis will be used. In most cases this hypothesis 
is stated as follows : During bending, plane cross sections which are 
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normal to the axis before deformation remain plane after deformation, and 
normal to the deformed axis. Usually, this hypothesis is combined with 
the assumption, although not always stated, that strains within the cross 
sections can be neglected. This assumption will be used in the present 
study also. (This is similar to the case of plate and shells where the 
analogous Love-Kirchoff hypothesis is used. ) This hypothesis leads to 
the following results: 


and 



R 


x,, e + W + y„ e' + z_ e' + 
u x 0 y 0 z 


’’(“o * y 0 



(b- 18 ) 


(b-19) 


The last term in Equation (B-19) represents small normal displacement 
which, as pointed out by Novozhilov (Ref. 12, p. 213), is a generaliza- 
tion of the warping function of St.-Venant torsion. This function con- 
tains only quadratic and higher degree terms in y^ and Zq and is 
assumed to be small, compared with typical cross-sectional dimensions of 
the rod. 

Substitution of the expression (B-6) for W into Equation (B-19), 
then differentiating (B-17) and using Equation (B-17), implies: 


G 

x 


= R 


(l+u )e+v e+w e 
,x X , x y ,x z 

y (-* e ' + t e ' ) + z_(-k e ' - t e * ) 
J 0 V y x z' 0 V z x y ' 


+ <P e’ + <p(k e' + k ) 
,x x y y z z 


(B-20a) 
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(B-20b) 


G 


= 5 

■ e’ + 

q> 

e* 

>y 

y 

»y 

X 

- R 

A, 

= e' + 




z 


X 


(B-20c) 


By definition, e ' is a unit vector in the direction of E . 

x x 

aw 

Defining as the strain of the elastic axis, and then using the 

first of Equations (B-9), implies: 


E ® (l + u )e + v e+w e 
x ,x x ,x y ,x z 


(1 + e) 


e' 

XX X 


(B-21) 


From Equation (B-2l), using the phthagorian rule and neglecting terms of 

2 

order £ compared to unity, implies : 


xx 


, X / & C. \ 

+ 2 ( \* + * J 


2 

»x 


(B-22 ) 


Substitution of Equations (B-20) into the expressions (A-7) of Appendix 
A for the strain components, and making use of Equation (B-22 ), implies: 


xx 


xy 


xz 


e - y K - z_ K + ffi 

xx *0 y 0 z ,x 


i (*. y - ^0 T + * V 


(B-2 3 ) 


? to . + y 0 T + ¥ % ) 


In deriving the expressions (B-23) use vas made of the fact that 

the quantities y Q * y , z Q * z , y Q t, z q t, <p x , «p * y , ? * 8 , «P >y , 

2 

and «P are of magnitude of strain and therefore less than £ . 


B.3 Force and Moment Resultants 

In calculating the stresses in the rod, use is made of the con- 
stitutive relations of the material from which the rod is made. In the 
present case it is assumed that the rod is made of isotropic Hookean 
material which is homogeneous for every cross section. On the other hand, 
in the present study the assumption that o = = 0, commonly used 

for slender rods, is made. However, according to the Bernoulli-Buler 
hypothesis, £ and £ should also be zero, and the vanishing of 

%r j ““ 

c r , c , e and £ simultaneously, is inconsistent with Hooke’s 
yy zz yy zz 

Law. This inconsistency, which iB inherent in the Bernoulli- Euler hypo- 
thesis, although not always stated is explained in the literature in 
different ways. One of the explanations is that one is dealing with a 
material having a special type of ort hot ropy. 

The constitutive relations for an orthotropic material, as given 
by Lekhnitskii (Kef. lb, Eq. (3-7)), are: 
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It is assumed, that in the present case: 


v 2l = v 3l = V 12 “ v 13 " 0; h " e 3 - " ; °23 58 G l3 " 2G * (B_2Ub) 

By using Equation (B-24b) together with the relations (B-24a), the incon- 
sistency that was mentioned earlier, disappears. 

Another inconsistency which is inherent in the Bemoulli-Euler 
hypothesis, concerns the shearing stresses. If torsion is neglected, 
the assumption that plane cross sections, before deformation, remain 

a# 

plane after phe deformation, means that the shearing strains £ and £ 

xz yz 

are zero. (These components are due to contributions other than torsion. 

The contributions due to torsion appear in Equation (B-23). ) This means, 
according to Hooke's Law, that the shearing stresses due to contributions 
other than torsion disappear. However, these stresses, and or 

more accurately their resultants — the shearing forces — do not vanish at 
all. Furthermore, they play an important role in the equilibrium calcula- 
tions . Sometimes this inconsistency is explained by taking Gt^ “ -* 

00 in Equations (B-24a). In the present case, where shearing strains 
due to torsion are also present, this assumption will cause scxne 
problems* Therefore, a better explanation is that in the case of slender 
rods the shearing strains are very small, so that they do not violate the 
hypothesis; however, their integral over the cross section should be taken 
into account, implying that the shearing forces cannot be neglected. 

Using Equations (B-24a), (B-24b) together with the strain relations 
as given by Equation (B-23), implies: 

*The two in the expressions for G^ and is needed because shearing 

strains in Ref. l4 are defined without the factor l/2 which appears in 
Eq. (A- 7) of the present study. 
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c = E E 
xx xx 


E(e - k - z n k + 9 ) 

v xx J 0 y 0 y ,x' 


T = T + 2G e 
xy xy xy 

= T + G(CP - J. T + 9 K ) 
xy ,y 0 y 


T «■ T + 2G E 
xz xz xz 


T + G(cp „ + y n T + <p K„) 

XZ | Z v z 


(B-24) 


Detailed expressions for and are not required, as will be shown 

later. 


The force which acts on the unit area of the cross section of the 
deformed rod, is : 


t 53 c G+t G+t G 
xx x xy y xz z 


(B-25) 


Using Equations (B-20), (B-2l), combined with (B-25), one obtains: 


[o(l+E -y K-zK+<p)+T Cp + T 9 ]«' 

xx xx J 0 y 0 z ,x xy ,y xz ,z x 


+ O (-z rt T + 9 K ) + V ]e' 
xx 0 y xy y 


+ [<r xx (y O T + * O + T -J 


xz z 


(B-26) 


Assuming that the stresses are of the same order and neglecting terns 

2 * 
of order e compared to unity, implies: 

See coament on page 57. 
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(B-2 7 ) 


t = a e’ + t e’ + t e' . 
xx x xy y xz z 

The resultant force, F, which acts on the cross section, is 
obtained by integration: 

F = // t dy_ dz n = T e* + V e' + V e' . (B-28) 

JJ ^0 0 x y y z z 

A 

Before proceeding, the warping function, which until now was treated 
in a general manner, has to be considered. One possibility it to treat it 
in an exact fashion as done by Wempner (Ref. 11, Chapter 8). This procedure, 
however, complicates the derivation considerably. Instead, when dealing 
with a slender rod, it is possible to introduce an assumption analogous to 
the one used for the case of the St. Venant torsion, whereby q> can be 
written as: 


* “ T ^ (x o»y 0 ' z o } 


(B-29) 


where <p is still a function of because it is conceivable that the 
cross section changes along the span with Xq. 

According to Equations (B-28), (b-2^), (r- 27) and (B-29): 


T 


rr ^ 

// E(e - y_ k — z_ k + t <p + t ® )dy_ dz„ 
JJ ' xx *0 y 0 z ,x T ,x ^0 0 

A 


(B-50) 


Restricting the derivation to symmetric cross sections (at least about one 
axis of symmetry) yields: 
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J7 5 ** 0 dz 0 

A 


II v % ^0 

A 


» o 


and Equation (B-30) becomes: 


T 


E A(e 


xx 


- y. 


oc y 


oc 


K 

z 


) 


(B-51) 


(B-32) 


while: 


II y o % to o ■ *.,= A > 

A 

(b-33) 

II *0 % *‘0 * Z oc A • 

A 


The point (y ,z ) is the point of intersection of the tension 
^ OC oc 

axis and the cross section. Equation (b- 32 ) together with the first of 
Equations (b- 24) implies: 

cr = S + (y • y.)E k + (z - z„)E k 

- - xx A J oc 0 y oc 0 z 

+ E t cp+ETCp . (B-34) 


After examining the force resultants on a cross section of the de- 
formed rod, the moment resultants about the point (y Q = z Q = 0), M, will 
he considered, where: 


M 



A 


d X t dA 


(B-55) 
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According to the Bemoulli-Euler hypothesis : 


3 “ y 0 ®y + Z 0 + T 9 *x * 


Substitution of Equations (B-27), (b- 24) and (B-36) into (b- 35) 
implies : 


(B-36) 


M » M e' + M e' + M e' 
xx y y z z 


(B-37) 


where : 


M =» / / [y_ t - z r ]dy_ dz_ 

X JJ xz 0 xy^O 0 

A 


M y ' JJ *0 - T "> T ra ldy o % 

A 


M 


JJ [ -V »b * T ? T w 1 % ^ 


(B-38) 


Substitution of Equation (B-24) into the first of Equations (b- 38) 
implies : 


M 


JJ (y o - *o \y )dy o te o * G S T JJ *0 5 % ^0 

A A 

; ° v t ir *o ? % % 

A 

+ G T iF [y 0 + Z 0 + y 0 - Z 0 *,y ldy 0 dz 0 * 


(B-39) 
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The first integral in Equation (B-39 ) Is the torque which is produced 

by the shearing forces V ,V around the point z ■ y =* 0. This point 

y z oc oc 

is the shear center of the cross section; thus, by definition this integral 
becomes zero. The last Integral is the torsional stiffness, J, of the 
cross section known from St. Venant torsion. Thus, expression (b- 59) 
simplifies to: 


M x - G J T + G T j/ y 0 ? dy 0 dz 0 

- 0 T V JJ z o ? % dz o • (B - to) 

A 

Substitution of Equations (B-34) and (B-24) into the second and 
third of expressions (b- 38) implies: 


M y * - E «y - S I33 \ + T 2oc 4 t >x E JJ z 0 V dy Q dz Q 

A 

+ T E JJ z 0 dz o - T JJ 5 \z % dz 0 

A A 

- t2<5 JJ % dz o - t2 g JJ *0 5 % % 

A A 


• T % g JJ ? % dz 0 ! 

A 


(B-4la) 
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rr 2 

hi - jj (y oe - V % dz o 

A 

hi - JJ <z oc - z o )2 % iZ 0 

A 

rr 

^3 " Jj (y oc - y o )(z oc - z 0 )dy o dz o 

A 

The underlined terms in Equations (B-^l) become zero in the case of a 
symmetric cross section, which is the type of cross section being considered. 

In the case of slender rods with closed cross sections, the influence 
of warping is usually neglected and in this case, expressions (B-4 o) and 
(B-^l) become: 
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M * G J T 
X 

M y * - E *23 V - E bi “z * T V 

E \z V + E ^3 "z - T y oc 

Comment : In Equation (B-26) a contribution of the axial stress 

to the shearing forces acting on the cross section of the blade exists. 
For beams where the torsional stiffness is very small compared to the 
bending stiffness (like the case of beams with thin open cross sections) 
this contribution, sometimes called the trapeze effect, can cause con- 
siderable influence of axial forces on the torsional rigidity of those 
beams (see, for example. Ref. 16). However, rotor blades, which are the 
subject matter of this study are made of either closed or solid cross 
sections, where the above mentioned effect can be neglected, as pointed 
out by Goodier (Ref. 16, p. 386, second column, line 18 from the top). 
Therefore, the assumptions leading from Equation (B-26) to (B-27) seems 
to be appropriate for the present study. 


(B-U3) 
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Fig. B3. Euler Angles When the Order of Rotation is 9 9 8 

z y x 
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APPENDIX C 


EQUATIONS OF EQUILIBRIUM OF A DEFORMED ROD 

In Appendix B the expressions for the resultant forces and moments 
which act on a cross section of a deformed slender rod, which was initially 
straight, were calculated. Furthermore, it was assumed that the rod is 
subjected to a distributed force, p, per unit length of its undeformed 
axis. This load p includes body forces, surface tractions and inertial 
loading. There is also a moment q per unit length of the undeformed 
axis of the rod. This includes body couples, moments of surface tractions 
and moments of inertial loading. The loads p and q are assumed to be 
continuous and also having continuous derivatives. 

Figure C-l shows a segment of the deformed rod. From equilibrium 
of forces the following equation is obtained: 

F + p - 0 . (C-l) 

From the equilibrium of moments about the point P, letting dx^ -» 0, 
the following equation is obtained: 

M + q + e 1 X F = 0 . (C-2 ) 

y X X 

Hie load p and the couple q are described by their components : 
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p 


p e' + p e' + p e' 
x x y y z z 


q 





(C-3) 


Substitution of the first of Equations (C-3) into Equation (C-l), 
together with Equations (B-17) and (b- 28), implies: 


T 

K V — 

K V + 

p 

= o 

(e' 

direction) 


t x 

y y 

z z 

X 


X 

j 

V + 

K T - 

T V + 

p„ 

= 0 

(e* 

direction) 

| . (C-4) 

' 

y,x 

y 

z 

y 


y 

V + 

z,x 

K T + 
z 

T V + 

y 

p z 

= o 


direction) 



Substitution of the second of equations (C-3) into (C-2 ), together 
with Equations (B-17), (b- 28) and (B-37) implies: 


M — k M - k M + q 

x,x y y z z 

= 0 

% 

direction) . 


M + * M - TM +q-V 

y>x y x z t z 

= 0 

% 

direction) 

. (C-5) 

M +K M + tM + q_ + V 

z,x z X y * y 

= 0 


direction) 



Equations (C-4) and (C-5 ) are exact, and contain no approximations. 
The procedure of solution is as follows: Expressions for and 

are obtained from the second and third of Equations (c-5). These are 
subsequently substituted in Equations (C-4). Following this procedure, 
the expressions for the shearing forces become: 
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-(M + k M+TM+q_) , 

' z,x z x y ^z ' 


M +KM -TM+q_ . 

“ z My 


y,x y x 


( 0 - 6 ) 


Substitution of expressions (C-6) into Equations (c-4) implies: 


T +k(m + t m + q_ 
,x y' z,x y 


) 


- K JK, „ - T + Qy) + P„ - 0 


z y,x 


(c-7a) 


-(M + k M + t M + q_) +K T 

z,x z x y ^z ,x y 

-t(m +kM-tm+o)+p 
y,x y x z > y 


(C-7b) 


(M +kM-tm + a ) + k T 

• y,x y x z t ,x z 

- t(m + k M + TM + q_) + p 
z,x z x y y. z 


(C-7c) 


Equations (C-7) represent the equilibrium of forces in the direc- 
tions e', e ' , and e', respectively. There is still the moment equa- 
x y Z 

tion, in the e ' direction that must be satisfied, which is: 


M -KM 
x,x y y 


K M + 
z z 


V 


0 


(C-8) 


Equations (C-7) and (C-8) are accurate and contain no approxima- 
tions. These four equations (three of (C-7) and one of (C— 8 ) must be 
solved in order to investigate the problem of the deformed rod. In order 
to obtain a solution, it is necessary to express the moments in terms of 
the derivatives of the displacements and the rotation of the cross 
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section about the elastic axis. This reduces the problem to one contain- 
ing four equations and four unknowns. 

In trying to simplify the equations, use can be made of the order- 
ing scheme. Performing the differentiation in Equations (C-7a - c) 
implies : 


T + * M -KM +t(k M+k m) 
,x y z,x z y,x y y z z 


+ q_-K q+p = o 

y T5 z y x 


(C-9a) 


-M - (k + tk )m - k M - t M 
z,xx z,x y x z x,x ,x y 


- 2t M +T 2 M+KT-q_ - T q_ + P =0 (C-9b) 

y»x z y H z,x T y 


M + ( K - t k )m + k M -t m - 2t M 
y,xx y,x z x y x,x ,x z z,x 


- T M + K T + 
X z 


V,* ~ T **« + P * 


0 


(c-9c) 


M -KM-«M+q_ = 0 

x,x y y z z tc 


(C-9d) 


The underlined terms in Equations (C-9b) and (c-9c) can be - _ . 


neglected according to the ordering scheme. As an example, consider 

2 

Equation (C-9b). The underlined terms t m can be neglected, com- 

z 

pared to the term -M which also appears in the equation. As a 

Z ; XX 

clarification, recall that according to the ordering scheme one can 
write: 


M - M (1 - cp 2 ) . 
z z 
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If both sides of the last equality are differentiated twice with respect 
to x> it implies: 


M 


z.xx 


M - M 

Z,XX Zj 


XX 


«p 2 - hM w 


z,x 


- 2M (f 
z , 


W ) 


XX 


The underlined terms in the last equality are negligible compared 

2 2 

to M . The term -2M «p is of the order of magnitude of t m 

Z f XX Z > X z 

and so neglection of the underlined terms in Equations (C— 9 ) is justi- 
fied. It is clear that the last argument is correct only if deflections 
and force and moment resultants are changing gradually along the span 
and do not have very high gradients. 

Substitution of the expression for M from Equation (c-9d) 

XjX 

into Equations (C-9b,c), and using the ordering scheme, Implies the fol- 
lowing set of equations: 


T +*M -KM +t(ic m+k m)+k 
,x y z,x z y,x y y z z y 


- k q + p =0 
z y *x 


(C-lOa) 


-M - (k + t k )m - (t + k k )m - 2t M 
z,xx z,x y x , x y z y y,x 


+ K T + K 

y z 


q - t q - q + p = 0 

tc y z,x y 


(C-lOb) 


M + (k . . _ T k )m - (r - k k )m - 2t m 

y,xx y,x z x ,x y z z z,x 

+ k T-k q - Tq+n + p = 0 

z y tc y y,x z 


M — k M — k M+q_ = 0 

x,x y y z z ^ 


(C-lOc ) 
(C-lOd) 
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Sometimes it is more convenient to write the equations of equi- 
librium in the undeformed directions (e ,e ,e ), instead of the 

x y z 

directions after the deformation ) as was done in the previous 

x y z 

part of this Appendix. The general relation between the two systems is 
given in the form: 


e* = e +S. 0 e + S, , e 
x x 12 y 13 z 


e' = S 0 . e + e + S„, e 
y 21 x y 23 z 


e' = S,.. e + S-_ e + e 
z 31 x 32 y z 


(c-ll) 


Two such transformations, belonging to the class of transformations 

2 

containing the inherent assumption that quantities of order of £ are 
negligible compared to unity, are presented in Appendix B as Equations 
(B-13) and (B-lU). The components are of order, e, or less. 

The resultant elastic force on a cross section of the rod is given 
in an analogous manner to Equation (B-28), by the expression: 


■ • F -■ T e + V e + V e- . ■ — . .. . . (C-12 ) 

x y y z z 

and the resultant elastic moment, which acts on a cross section of the 
rod is given in an analogous manner to Equation (b- 37), by the expression: 


M = Me+Me+Me 
xx y y z z 


(c-13 J 


Equations (C-12 ), (C-13)> (B-28) and (B-37) together with Equation 
(C-ll) imply: 
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1 + S 21 V y + S 31 V s 


S 12 T + V y + S 32 V z 


(0-14) 


and 


S 13 T + S 23 V y + V z 


M = M + S_, M + S x , M 
x x 21 y 31 z \ 


M = S,~ M + M + S,„ M 
y 12 x y 32 z 


M “ S,M + S„ , M + M 
z 13 X 23 y 2 


(C-15) 


The distributed force p and distributed moment q per unit length 
are given by (compare to Equation (C-3)): 


P " P x e x + P y e y + p 2 e 2 


q x e x + q y e y + e * 


(c-16) 


From the equilibrium of forces (Equation (C-l)) the following equa- 
tions are obtained: 


T + p 
,x x 


V + p 
y,x y 


v + p 

z,x z 


(C-17) 


and the equilibrium of moments (Equation (c-2 )) Implies: 
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0 


M ' M M m 


M - S-, 

x,x 15 

V + V + 

y 12 z 


» 0 

i 

1 


M + S., 

y»x 15 

+ 

i 

IH 

V 

- 0 

> 

(c- 18 ) 

M - S 

z,x 12 

T + V + 

y 

V 

= 0 




From Equations (C-14), neglecting terms of order e compared to 
unity, one obtains: 


5 - I+ (s 2i - V23 )7 y + (S 31 - S 21 S 32 )7 * * (c - 19) 


The second and third expressions of Equations (C-l8), together 

2 

with (C-19), yields, after neglecting terms of the order e compared 
to unity: 


V - -M + - S„-S,„)M 

y z,x 12' 51 21 52 y,x 


+ S 12 T - \ + S 12 (S 51 “ S 21 S 52 } V ' (C ' 20a) 


and 


V - M - S. X (S_. - S,-S q ,)m + S,, T 
z y,x 15 21 51 25' z, x 15 


" S 15 (S 21‘ S 51 S 25^ q z + V * 


(C-20b) 


Substitution of Equations (c-19) and (C-20) into (C-17) and the 
first expression of Equations (C-l8), yields: 


[T ‘ (S 21 " S 51 S 25^ M z,x + (S 51 “ S 21 S 52^y,x ] ,x 


- I(s 5L“%1 8 32 , V 1 „ + 5 ? x = 0 ’ (c - 21a) 
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(C-21b) 


[ - 5 z,x + S 12 (S 31 - S 21 S 32 y,x^,x + (s !2 *)„ 


- Vx + [S 12 (S 31 ' S 21 S 32^y',x + ? y 


0 


[ *V,x ' S 13 (S 21 " S 3l S 23 )M z,x ] ,x + (S 13 T) ,x 

+ V,x “ [S 13 (S 21 - S 3l S 23 } ^ ] ,x + *z “ 0 


(C-21c) 


M + S_, M + S „ M + S„ a + S. x a + q = 0 . (C-21d) 

x,x 13 z,x 12 y,x 12 13 ^ tc 


According to the assumption that quantities of the order of e are 
negligible compared to unity, and the orthanormality of the transformation 
represented by Equation (C-ll), one has: 


21 " 

S 3l S 23 

* 

- S !2 

31 ' 

S 21 S 32 


- S 13 

+ 

s, 

s 

-s„ 

12 

31 32 


21 

4* 

s s 

rw 

_s. 

13 

21 23 


31 

23 “ 

S 12 S 31 


ft 

W 

1 

32 “ 

S 13 S 21 

fW 

- S 23 

21 + 

S l3 s 23 


' S 12 

31 + 

S 12 S 32 


" S 13 




(C-22 ) 
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Making use of* Equations (C-2 3), Equations (C-21a-d) turn out 


to be : 


(T+S 12 5 «,*- S 13 5 y,* ) ,* 




(G2 3a) 


<A,x - S 12 S 13 \.x\x + (S 12 T) ,x - 
- <S 12 S 15 V,x + *y = ° 


(C23b) 


(M + S,_S, , M ) + (S_, T) + a . 

' y,x 12 13 z»x ,x 13 ,x T»- 


* ‘VlS^.x* P * 


0 


(C23c) 


m + s,, m + s _ m + s -,o <1, + s ,i <l + %■ 

x,x 13 z,x 12 y,x 12 T 15 ^ * 


0 . (C-23d) 


Substitution of Equations (C-15) into Equations (C-23), using 
Equation (C-24d) also, to substitute for M x , and neglecting terms 


of order e compared to unity, yields: 


! T - S 21 v* + S 51 M y,x - S l5 ( V,x m x* W.x H y 

♦ tS I2 (S 13 ) ,x- S 15 (S W ) ,x lM x 5 ,x + (S 12 %>,x- (S 13 %\x + ? X - 0 (C - 2 ‘“ ) 


■ (M I,X + V t(s 2 3 U ■ s 13 (s 21 ) ,x lM y " s 52 M y,x ! ,x 

+ < S 12 T) ,x - V* (S 15 ^x* *y ’ 0 (C ' 2lrt,) 
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M 

x, x 


+ (S 13 T) ,x - (S 12 \\x + «y,x + ** = ° ’ 


* t( V,x + [( V,x + W,x lM 


+ Sc + S 12 V + S 13 % 


( C-2 4c ) 


z 


0 


( C-2 1+d ) 




APPENDIX D 


DERIVATION OF THE EQUATIONS OF EQUILIBRIUM BY THE 
USE OF THE PRINCIPLE OF VIRTUAL WORK 

D.l Principle of Virtual Work Applied to a Rod 

Usually, in a structural system subject to loads, internal forces 
develop between the various components of the system as a result of the 
external loads. If the work of the internal forces is denoted by Wj, 
and that of the external forces W £ , then the principle of virtual work 
can be simply expressed as: 

sw I - - ew £ (D-1 ) 

where 9W^ and SW E are the work done during a virtual displacement by 
the internal and external forces, respectively. 

In the case of an elastic system. 


U 



(D-2 ) 


where U is the elastic energy in the system. Therefore, Equation (D-l) 
becomes : 


8U = 0W„ 

E 


(D-3) 


Using Equations (D-2 ) and (D-3) implies (for example, Eq. (9-125) 



of Wempner, Ref. ll), the following expression for a continuous "body: 


Wj - 


Volume of the 
body 


[<r 8£ +a 8E + c 6e 

xx xx yy yy zz zz 


+ 2T 8E + 2T SE + 2T 8E ]dx~ dy_ dz_ 
xy xy xz xz yz yz u J 0 0 


yz yz 


For the case of a slender rod of length i and cross section 


A, within the framework of the Bernoulli-Euler assumptions 


a = ci =e =o, 
yy zz yz * 


Equation (d- 4) becomes: 



8E + 2t 8e + 2t 
xx xy xy xz 


se xz ld:, b % • 


The virtual displacement of the rod is given by a displacement 
8W (xq) of every point on the elastic axis, accompanied by rotation 
8© (xq ) of the triad e^,e^,e^ at every point of the deformed rod. The 
virtual rotation 8© can be described in this case as a vector because 
it is infinites imal - The rod is acted upon by a distributed force, p, 
per unit length of its undeformed axis, which includes body forces, 
surface tractions and inertia loading, and it acts at the elastic axis 

The fact that £y Z * 0 emerges from the assumption that strains within 
the cross sections of the rod are neglected. 


(d-4) 


(d-5) 
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as described in Appendix C. There is also a distributed moment, q, per 
unit length of the undeformed axis of the rod. It includes body couples, 
moments of surface tractions and moments of intertial loading. This 
moment also acts at the elastic axis of the rod. It is clear, therefore, 
that : 


8W 


_ r 


*0 


= l 


E 


(p • 6W + q • 5©)dx 


*0 


= 0 


0 


(D-6) 


Equations (D-3), (D-5) and (D-6) imply: 



// [a be + 2T 8£ + 2T 8£ ]dv„ dz 

JJ xx xx xy xy xz xz y 0 0 

A 


p »8W - q *89 


% 


0 


(D-7) 


According to Equations (B-23) of Appendix B: 


8£ 


XX 


6E 

xy 


be 

xz 


5£ - y. 5 k 

xx J 0 y 


- 2 z 0 5T 


l y 0 6T 



(D-8) 


In deriving Equation (D-8), use was made of the assumption that 
in the case of a slender rod with closed cross section, the influence of 
warping is negligible and, therefore, the virtual work is taken as that 
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performed during a rigid body motion of a slice of the rod, and thus, 8cp 

is assumed to be zero. The terms <p8k and <P8 k are also neglected as 

y z 

a consequence of assuming the warping to be negligible. 

Substitution of Equation (d- 8) into (D-7), use of Equations (B-38) 
from Appendix B, and neglecting the influence of warping, Implies : 

r *o ui 

J (T6£ xx + M z 8K y- V\ + M x 6T - q.89)dXo = 0 , (D-9 ) 

" X 0 = 0 

where 8£ , 5k, &k and &t are functions of 8W and 60. After 

xx y z 

substitution of the appropriate values for 8£ , &k , 8k and 6 t, the 

xx y z 

equilibrium equations and boundary conditions are obtained by performing an 
integration by parts of Equation ( D-9). 

In Chapter 2 of this study it was shown that the equilibrium 
equations can be obtained with respect to different directions. In the 
present case, the definition of &W and 80 determines the directions in 
which the equations will be valid. 



D.2 Equilibrium Equations in the Directions of the 

Deformed Rod Coordinates (e' ,e' ,e' ) 

x y > z ' 

In this case, the virtual displacement is chosen as: 

&W = 8u' e' + 8v' e ' + 6w' e' . (D-10) 

x y z 

The virtual rotation is given in the form: 
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'> a A. 

n e' + n e' + n e' 
xx y y z z 


(D-ll) 


where n ,n and n are the components of the virtual rotation. 
x y z 


Due to 6W and &© the triad e'.e'.e' is rotated to a new 

x ' y 7 z 


triad e£ , e^. , e^ , « iven : 


e" • e' + 6© x e’ 

XX X 


e' + n e' - n e’ 
x z y y z 


e' + 6© x e* 

y y 


A. '' f 

-n e' + e' + n e 
z x y x z 


(D-12 ) 


e' + &© x e' 
z z 


“ n e' 


A | Ss 1 

n e' + e 


y x x y 


It is clear that n = & <f> , however n and n axe determined by &W. 

x y z 

In order to find n y and let us consider an element cLXq of the 

deformed elastic axis, which is shown in Figure D-l. Before the virtual 
displacement, the element is in position AB, described by dXg e^. 
After the virtual displacement, the element is in position A'B’, given 
by : 


A'B’ = (R + dxQ e^ + BW + 6W x dxQ) - (R + &W) 





(D-13) 


Substitution of expression (D-10) into (D-13), and performing the differ' 
entiation, while using Equation (B-17) of Appendix B for the derivatives 
of the unit vectors, implies: 
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(D-14) 


A'B' = [ (l + 5u' -* 6v' - k 8w')e' 

,x y z x 

+ (&v* + K 5u' - T 8 w' )e' 

y y 

+ ( 8w' x + 8u' + t 8v' je^ldxQ 

Recalling that the virtual displacement 5W is as small as desired, then 
from Equation (D-l4 ), one obtains after neglecting products of virtual 
terms : 


^ M 


e 


x 


e' + ( 8v' + * 6u' - t 8w')e' 

x ,x y y 


+ ( 8w' + k 8u' + t 5v')e' 

,x z z 


(D-15) 


Comparing Equation (D-15) with the first of Equations (D-12 ), 

the quantities n and n are determined. Thus, the rotation components 
y z 

are : 


n = 8 » 
x 


n = -( 5w' + k 8u* + t 6v' ) 

y ,x z 

n = ( 6v' + k 8u' - t 8 w' ) 

z ,x y 



Using definitions 


(B-15 ) from Appendix B, then: 


(I>l6) 
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K + B*C 

y y 


K + 8* 

z z 


e 


x.x 


e " 
x,x 


e" 

z 


(D-17) 


T + 8t 


e " 
y»x 


Differentiation of expressions (D-12 ), using Equations (B-17) 
of Appendix B, yields : 


e " 
x,x 


c e' + k e ' + (n * - n * )e 1 

yy zz yz z y x 


+ (n + Tn )e' + (-n +t n )e' 
z,x y y y,x z z 


(D-l8a) 


e" ■ -k e' + t e' - (n + n k )e' 
y,x y x z z,x x z x 


- (n t + n * )e' + (n - n * )e' • (D-l8b) 

x z y y x,x z * z 


Substitution of Equations (D-l8) and (D-12 ) into Equation 
(d- 17), and neglecting nonlinear terms in n ,n and n , implies: 

^ y * 


6* 


Bit 


n + k n + t n 
z,x z x y 


-n - * n + t n 
y,x y x z 


(D-19) 


B T 


n - k n. - * n 
x,x y y z z 


Substitution of expressions (d- 16) into (D-19) yields: 


Bk » 8 v* + (k 8u') - (t Bw' ) 

y >xx y fX ,x 


+ K 8^ - t Bw' - r K Bu' - T Bv' , (D-20a) 

z ,x z 
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9 


(D-20b) 


8 is = Bw' + (k &u’ ) + (t 6 v' ) 

z ,xx z,x ,x 


- K 8 $ + T 8 v’ + T K 8 u' - T gw' 

y ,x y 


gT = g 1 ® + K gw' + K K 6u' + K T 6v' 

*X y ,x y z y 


K 6 v' - K K 5 u' 4 - K T 8 w' 

z ,x y z z 


(D-20c) 


From Equation (D-l^ ) it is clear that (neglecting products of components 
of the virtual displacement ) : 


ge = & u' - k 6 v' - is gw' . 
xx ,x y z 


(D-21) 


The load p and couple q are given vith respect to the de- 
formed system, in the form of Equations (l) and (2), as: 


p e' + p e' + p e' 
x x y y z z 


* 3 V*y + 


(D-22) 


Substitution of expressions (D-10), (d-11), (d-16), (d- 20), 
(D-2l), and (D-22 ) into Equation (l>-9) implies: 


r o 

/ (T 6u' - T K 8v' - T is Bw' + M 8v* + M (is 8u') 

J ,x y z z ,xx z y 

•*0 

- M (t gw' ) + M is 8 ♦ - M T gw' - M t x 8u' - M t 2 8v' 

z ,x z z z ,x z z z 

- M Bw' - M (is Bu' ) - M (t &▼•) + M k 6* - MJ 8v' 

y , xx y' z fX y' y y y 
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- M T K 

8u' 

+ M 

T 2 6w' + 

M 6 ♦ + M K 

8vr ' 

+ M K K 

8n' 

y 

y 

y 


x ,x x y 


x y z 

+ M 

K T 

8v' 

- M 

K 8 V* - M K K 8u* + 

M K 

T 8w' 


X 

y 


X 

z ,x 

x y z 

x : 

s 


- p x 

8u* 

- p 8v' 

y 

- p 8 v* 
z 

- + q,**;,* 

a k 8u* 

y z 

+ q^T 5v* 


- q^ K^8u* + \ r &' r ' ) dx o “ 0 • (D-23) 


Integration by parts of Equation (D-23) gives the following 
variational expression: 



(-1^ 5u' - Rg 8v’ - R 5 8w* - R^ 6 ♦ 


+ 

[b Su 1 1 0 + [a 

-^*0 J *b“° L 



L1 Jxq^O L 

-v° 


r -i x °“ i 

II 

>? 

r 


+ 

R. 8v ' + 

J *Xj r\ 

^ * - j Xq *0 

L E 6 *^0 * 0 

(D-24) 


where the various expressions in (d- 24) are: 


R. =* T +k(M +tM+q)-k(M -TM+q_)+p , (D-25a) 

1 ,x y z,x y ^z z' y,x z ^y *x * 

It, = -(M + K M + tM + a) + k T 

^ z,x z x y ^z ,x y 

- t (M + k M - TM + q. ) + p , (D-25b) 

y,x y x z ^y *y ’ 

R, = (M +KM-TM+q) +KT 
3 y,x y x z y ,x z 

- * "z M X * ™y + + P« > (D-^c) 
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and 


r = M -km 
4 x,x y y 


k M + 
z z 


(D-25d) 


B, - T+k M - k m 
1 y z z y 


EU * - (M +2tM+k M+q) 

c z,x y z x H z 


B, = M +KM-2TM + 
3 y,x y x z 


B,. = M 




B,- = M 


B c - -M 
6 y 


(D-26) 


With the virtual displacement being arbitrary, the equilibrium 
equations turn out to be: 


= 0, Rg = 0, R^ = 0, - 0 • ' - (D-27) 

Substitution of Equations (D-25) into (D-27) yields exactly the same 
equations of equilibrium (C-7,8), which were obtained in Appendix C. 

The boundary conditions for this problem can be obtained directly 
from Equation (d- 24), subject to the assumption that the boundary condi- 
tions are not varied during the loading. 
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B 1 = 0 or 
Bg = 0 or 
B^ = 0 or 
B^ = 0 or 
B._ = 0 or 
B^ *■ 0 or 


u = 0 
v = 0 



(D-28) 


The boundary conditions of free edge and rigid clamping, as 
stated by Equations (9) and (10) of Chapter 2, are in agreement with 
Equations (D-28). There is only one item which should be noted. In 
Equation (9) the boundary condition of u = 0 at x Q = 0 is neglected. 
This is due to the fact that T is used in the equilibrium equations, 
instead of u. As a result of this difference the first equation con- 
tains the first derivative of T instead of a second derivative of u. If 
expression (B-32 ) of Appendix B is inserted in the equilibrium equations 
instead of T, using Equation (B-22 ) for £ , the boundary condition, 

XX 

u = 0 at Xq » 0, is needed when the unknown T is replaced by u. 


D.3 Equilibrium Equations in the Directions of the 
Undeformed Coordinates of the Rod (e ,e ,e ) 

x > y> z 

As pointed out in Appendix C (Equation (C-ll)), in general: 


8k 



(D-29) 


e' = e + S, 0 e + S, , e 
x x 12 y 13 z 


^ . _ A - _ A 

e' = So, e + e + S ox e 
y 21 x y 23 z 


A , _ A « A A 

e' = S,, e + S,o e + e 
z 31 x 32 y z 


Because of the orthonormality conditions. Equations (d- 29 ) imply: 


= e' + s p , e' + S,_ e' 
x 21 y 31 z 


_ . A, A 

S„o e' + e* + S, 0 e 


12 x 


32 z 


(d*30) 


S, x e' + e' + e* 
13 x 23 y z 


The curvatures and twist, using definitions (B-15) of Appendix B are: 


k = 


t = 


A, 

e' • 

e* = 

(S-o) 

x,x 

y 

12 ,x 

a, 

e' 

e’ = 

(S.J 

x,x 

z 

13 ,x 

A. 

e “• 

y>x 

A| 

e = 
z 

(S 23 ) ,x 


(D-31) 


In this particular case, the virtual displacement is chosen as 
(compared to (D-10)): 

8W = Bu e + &v e + 6w e . (d- 32) 

x y z 

Using Equations (D-13), (D-30) and (D-32 ), implies (compared with 

(D-l *)): 
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t-n * — 


A’B 


e' + Bu e + Bv e + 8w e 
x ,x x ,xY ,X Z 


■ (l + 8 “,x + S 12 8t ,x + S 13 6w ,x^ 

+ <S 21 8u ,x + 8 \x + S 23 8w , x )e ~y 
+ (S 31 8 ",x + S 32 6v ,x + 6w ,x )e l 


Neglecting virtual terms compared to unity in Equation (D-33 ), implies 
(similar to (D-15)): 

e “x * e ~x + (S 21 8u ,x + 8v ,x + S 23 Bw ,x^y 

+ (S 31 6u ,x + S 32 6t ,x + 8w ,x ) “z • 


The virtual rotation around the elastic axis is 8 0 e', thus, 

x 7 

according to Equation (D-12 ), one obtains: 


n =6-0 
x 


-y = - (S 31 8u ,x + S 32 8v ,x + 6w ) x > 


m = S„_ Bu + 8v + S oz 6w 

z 21 ,x t x 23 ,x 


Using (D-19)> together with (D-3l) and (D-35 )> implies: 


8K y =6l, ,xx* (S 21 8 U ,x ) ,x - < S 23 ) ,x S 31 8u ,x + 


' [(S 23 ) ,x + (S 2l\x S 31 ] 8v ,x + 1(s l3>,x + (V.x^ 1 


(d-33) 


(d-34) 


(d-35) 


5 0 (D-36a) 
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(D-36b) 


^z = Bw ,xx + (S 31 8u ,x\x + (s 23\x s 21 Bu ,x + (S 32 &V ,x } ,x 

* (S 21 J ,x S 51 1&r ,x- V Kj^x S 25 i8 * 


W ‘ S \x + [( V,x (S 5l‘ S 21 S 52 ) ' ‘ S 15>,x (S 21 - “ajV 1 6 “,x 


' {S nKx 8v ,x* (S 2Z‘,X 6w ,x 


. (d-36c) 


From Equation (D-33), one obtains: 


8£ = 6u + S 10 6v + S,,6w 

xx ,x 12 ,x 13 ,x 


(D-37) 


The virtual rotation components are given by Equations (D-35) with 

respect to the rotated triad (e',e' ,e' ). Using Equation (D-29), implies: 

x y z 


69 ■ [5. - (s 21 - S 2J S 51 )6 V+ (S 31 - S 2lS32 )6V jX ]e x 

+ [- 6w jX - ( s }1 - s 52 s 21 ) 8u jIC + 8»le y 

- [ 8v ,x + (s ai - 6u ,x + s « • <^ 58) 


With Equation (D-29) and the assumption that terms of order £ are 
neglected compared to unity: 


S 21 " S 23 S 31 = “ S 12 * 


S 31 " S 32 S 2l = " S 13 * 


(D-39) 
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The load p and moment q are given with respect to the undeformed 
triad (® x >®y>® z )> as Equation (C-l6) of Appendix C: 

— ** A A A 

p“pe+pe+pe , 

x x y y ^2 z 

(d-4o) 

i - Sc «x + V ®y + ^ ® z * 

Substitution of Equations (D-32 ), (d- 36 ), (D-37), (d- 38 ) and (l>- 4 o) 
into (D“9)» and using (D-39), implies, after integration by parts: 



" 0 (d-4i) 

where : 

h " ' T - S 21 M z,x + S 31 M y,x - (S 2}>,X S 21 M y 

+ M x [(S lj\x S 12 ' (S 12 ) ,X S 15^,x 

+ p - (S_ a ) + (s. 0 a ) 

y x ' 13 V ,x 1 2 ^z'jx 
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(S 12 T) ,x + % * <S 13 ^.x - V 


, (D-42b) 


£ M y,x + (S I z\x V [(S 23\x + (S 21^,x S 31 ]M z" S 23 M z,x^,x 
+ <S 13 T) ,x + \ + V,* ’ (S 12 *x\x 


«*,* - « s 15 >x* 


+ %c + S 12 V + S 13 "z 


B, 


% 



■ T - S 21 H *,x + S 31 M y,x- (S 23 ) ,x S 21 V (s 23 ) ,x Sx\ 
* M x [ (S 13\x S 12 ■ <S 12 \x S 13 ] ' S X3 V + S 3Z ** 

H z,x~ (S X}\x M x~ t(s 23 ) ,x ' f (a zi\x S 31 lM y 
+ s 32 M y,x + (S 12 T) ,x -%+ S U*x 

’ M y,x + (S 12 ) ,x M x- t(S 23 ) ,x' f (S 21 ) ,X S 31 lM z 


, (D-42c) 


, (D-42d) 


, (D-43a) 


, (D-43b) 


- S 23 M x,x + S 13 T+ V S 12 ^ 


, (D-U3c) 



> 


(D-43d) 


' S 3l M y + S 2l M z ' 


-S^, M + M 
32 y z 


-M + M 
y 23 z 


The virtual displacements are arbitrary, thus, the equations of 


equilibrium become : 


\ * 0; Bg - 0; R, - 0; R^ - 0 j 


(l>-44) 


while the boundary conditions (assuming that the boundary conditions are 
not changed during loading) are: 


B x - 0 


u “ 0 


v *» 0 


°0 


\ - 0 


w » 0 


4> B 0 


(D-45) 


b 5 =o 


U = 0 
,X 


b 6 "° 


v =0 
,x 


B 7 - 0 


w “0 
,x 


Hie equilibirum equations (d-44) are identical to (C-24) of Appendix 
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axe 


2 

C, within the framework of the assumption that terms of order e 
neglected compared to unity, which implies: 


-s„- M - S.,(S TO ) M 

21 z,x 13 32 ,x z 

a 

“ s 2i M Z,X * S 31 (S 23\x M z ' 

(D-U6a) 

s 3i M y ,* + s i2 (£ W,* M y 

** 

S 31 M y,x - S 2l( S 2 3 } ,x M y ' 

(D-U6b) 



[(S 23^,x " S l3^ S 21^,x 5 » 

(d-46c) 

■^ s 15^,x + ^ S 12^,x S 3Z-' 

IV 

[ ^ S 3l^,x + S 12^ S 3e\x 5 * 

(D-46d) 

■ [(s 12 ) ,x + V.xV 

xs 

[ ^ S 2l^,x + S l3 (S 23\x 1 

(D-k6e) 


From the boundary conditions (D-V?) it seems that due to the change 

of the coordinate system a new boundary condition appears. B,. . ** 0 or 

u =0. However, upon checking the sixth and seventh condition, it is 

clear that in the case of a clamped edge, v ® w a 0, or a free 

,x ,x 

edge, My B ■ 0, the condition B,. » 0 is satisfied automatically. 
Thus, the fifth condition is satisfied, by satisfying the other boundary 
conditions. “ 
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Fig. D1. Displacement of Element on the Elastic Axis During the Virtual Displacement 
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